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Abstract 

Relation between Renyi entropies and moments of the Wigner 
function, representing the quantum mechanical description of the M- 
particle semi-inclusive distribution at freeze-out, is investigated. It is 
shown that in the limit of infinite volume of the system, the classical 
and quantum descriptions are equivalent. Finite volume corrections 
are derived and shown to be small for systems encountered in rela¬ 
tivistic heavy ion collisions. 


1 Introduction 

It is widely recognized that information on the entropy of the systems pro¬ 
duced in high-energy collisions is very useful in understanding the physics 
of the process in question. This is particularly important for heavy ion col¬ 
lisions and search for quark-gluon plasma. It was suggested some time ago 
PP that measurements of coincidences between the events observed in high- 
energy experiments may provide an estimate of the Renyi entropies j2] of the 

*Address: Reymonta 4, 30-059 Krakow, Poland; e-mail:bialas@th.if.uj.edu.pi; 
^Radzikowskiego 152, Krakow, Poland 


1 



final state and thus also access information on its thermodynamical (Shan¬ 
non) entropy. Although the idea seems attractive, the argument of PP, being 
essentially of classical nature, can only be considered as a hrst step. The 
proper formulation must take into account the quantum-mechanical nature 
of the problem. In the present paper a quantum-mechanical formulation is 
developed and its consequences discussed. In particular, it is shown that the 
classical approach of PP is valid in the limit of very large size of the system in 
conhguration space. For systems of hnite size, quantum corrections are de¬ 
rived and shown to be relatively small for hnal states occuring in relativistic 
heavy ion collisions. 

The order I Renyi entropy, H{1), of a statistical system is dehned as 

= ( 1 ) 

where C{1) are coincidence probabilities of the states of the system, given by 

C{1) = Elfil' = Tr{pl (2) 

i 

The sum runs over all states of the system. Pi is the probability of a state i 
to occur and p is the density matrix of the system^. 

The attractive property of Renyi entropies is their relation to the standard 
(Shannon) entropy of the statistical system. It is easy to show that 

S = Tr[p\ogp] = \imH{l), (3) 

where S is the Shannon entropy. 

Moreover since, as is well known j2], for I > 1 

S >H{l)>H{l + l), (4) 

the Renyi entropies provide an exact lower limit for S, a quantity very im¬ 
portant for understanding the properties of the quark-gluon plasma 

The object of our investigation is an M-particle statistical system, i.e., a 
collection of M-particle final states which we define as those in which exactly 
M particles were observed in a given region of the momentum space. We 

^The second part of this equality is best seen in the representation where the density 
matrix is diagonal. 
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shall call them M-particle events (independently of how many particles were 
actually produced)^. 

On the classical level, M-particle hnal states can be described by the nor¬ 
malized M particle phase-space distribution K) with X = Xi, Zm-, 

A'= 

To obtain a quantum mechanical generalization of the phase-space distri¬ 
bution we follow the standard procedure, where the proper quantum descrip¬ 
tion of an M-particle final state is given by the density matrix p(pi,..., Pm; v'l) ■■■iV'm)- 
As shown by Wigner j7] the quantum-mechanical analogue of the classical 
phase-space distribution (called Wigner function), can be defined in terms of 
the density matrix as the Fourier transform 

W(X, K) = I dqe-‘'<^p(p, p'), (6) 

where K = [p + p')/^] q = p — p'. The quantum-mechanical description of 
multiparticle events is obtained by considering the Wigner function W{X, K) 
instead of the classical phase-space distribution hF^^“^(X, X). 

The main goal of this paper is to discuss the relation between the co¬ 
incidence probabilities C{1) defined in (j2I) and the moments of the Wigner 
function 


C{1) = J d^^X j d^^K[W{X,K)]K ( 6 ) 

Such a relation exists because both C{1) and C{1) are dehned in terms of the 
density matrix. 

The factor is introduced in Q to account for the factor 

(27r)“^^ in @ (which insures consistency between the proper normalization 
of the density matrix and of the Wigner function)^. This choice guarantees 
that, when hF(X, X) in (jHI) is replaced by the classical density distribution 

^This terminology is often used in experimental descriptions of multiparticle processes. 
For the momentum distributions, the proper technical terms are: exclusive distribution if 
all the particles are observed, and semi-inclusive distribution if besides a given number of 
observed particles there is an unspecified number of other particles. The latter should not 
be confused with inclusive M-particle distributions. 

^Even at the classical level, however, the phase-space distribution of particles produced 
in high-energy scattering is not a precisely defined quantity: one has to take into account 
that particles may be produced at different times. In the present paper, following ISIEI, 
we are considering the time-averaged distribution. 

Actually, these factors involve 27r?i. As is customary, we have put = 1. 
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{X, K), C{1) reproduces the correct classical limit of the coincidence 
probability, thus satisfying an important consistency constraint. We have 




:j 3 M 


\ 3 MTjrclass I 


l-l 


(7) 


In this formula /T) can be interpreted as probability distribu¬ 
tion and i^) as the probability for the system to occupy 

the elementary phase-space cell [271 around the (6M-dimensional) point 
{X,K). Thus one sees that, indeed, ([Tj) represents the classical formula for 
the coincidence probability. 

The interest in the relation between C{1) and C{1) stems from the obser¬ 
vation that the moments C{1) are experimentally accessible: we have shown 
recently jH] that, for a rather large class of models, C{1) dehned above can 
be estimated from the measured coincidence probability of the events 

with M particles [H El ITH] 




Ni 

N{N -l + l)/lV 


( 8 ) 


where Ni is the number of the observed 1-plets of identical events and N is 
the total number of events. N{N — l)...(iV — I + 1)//! is the total number of 
1-plets of events®. One sees that measurement of reduces to counting 

the coincidence between observed events. 

We show in this paper that C{1) (which can be measured by counting the 
number of identical events i), and C{1) (which determine Renyi entropies 
and thus open a window to the true entropy) are equal in the limit of inhnite 
size of the system. Also the hnite volume (quantum) corrections are studied 
and shown to fall with inverse square of the smallest (linear) size. 

In the next section a model for the Wigner function is introduced and the 
corresponding formulae for the moments C{1) are written down. In Section 3 
the results of jH] are summarized. In Section 4 the coincidence probabilities 
C{1) are analyzed in the same framework. The relation between C{1) and 
C{1) are discussed in Section 5. Our conclusions and outlook are given in the 
last section. 

®For 1=2 formula was first suggested, in a different context, by Ma El See also 

m 
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2 Moments of the Wigner function 

In terms of the Wigner function W{K^X)^ the momentum distribution is 
given by the integral 

w{K) = = J d^^XW{X,K). (9) 

To discuss C{1), we follow the argument of [S] and consider the time- 
averaged Wigner function of the rather general form, valid in a large variety 
of models 

W(X, K) = /■ G[XIL\w(K) (10) 

with X/L = {Xi — Xi)/Lx, {Zm — Zm)I Lz, K = Ki, Km- The function 
G satishes the normalization conditions 

J G{u)d^^u = 1^1 dXG{X/L) = {L^LyLz)^] 

J UiGiu)d^^u = 0 ^ < W, Yi, Zi >= W, % Zi; 

J{uifG{u)d^^u = 1 ^ >=T2,... (11) 

The hrst condition insures that w{K) is the correctly normalized (mul¬ 
tidimensional) momentum distribution, the second dehnes the central val¬ 
ues of the particle distribution in conhguration space and the third dehnes 
Lx, Ly, Lz as root mean square sizes of the distribution in conhguration space. 
Both sizes and central positions may depend on the particle momenta®. The 
form of function G is responsible for the shape of the multiparticle distribu¬ 
tion in conhguration space^. 

Using m we obtain from Q 

C(l) = (27r)“l'-‘> f [ -TTW|g(A7L)|' = 

= (2x9,)“<'-■) / 

®They may also be different for different kinds of particles. 

^The distribution in configuration space is given by f dKW{X,K) [J] and thus it is 
not identical to G(X/L). 
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with 


( 9 ,)«'('-!) = / d?'‘’u[G(u)]‘ 


(13) 


The constant gi depends on the shape of particle distribution in conhg- 
uration space. For example, we obtain for Gaussians 

and {gi)~^ = 2^/3 for a rectangular box. In the following we shall assume 
that G{u) is a Gaussian^ 


G{u) 


(27r)3M/2 


E m 

m = l ^-^a=x,y,z 


2/2 


(14) 


Introducing (HI into m we obtain 


C{1) 


/3M/2 J {L^LyL,)il-l)M • 


(15) 


3 Measurement of the moments of the Wigner 
function 

In this section we explain how one can estimate the moments of the Wigner 
function (ED by counting the coincidences (jHD between the measured events. 
The argument is a short summary of the results obtained in |H], where the 
relation between G{1) and G'^^^{1) was studied assuming the Wigner function® 
in the form given by eg. 

The major problem in the analysis of coincidences between the measured 
events is that these events are described by particle momenta which are 
continuous variables. Therefore, the dehnition (jSD is not directly applicable: 
a binning is necessary. Once events are discretized, the identical ones can 
be dehned as those which have the same population of the predefined bins. 

®This restriction can be avoided at the cost of some complications of the algebra. Since, 
however, the exact shape of the particle emission region is not well determined and since, 
moreover, o is not in obvious disagreement with the data from quantum interference, 
we shall stick to it. 

®In 13 the form llOII was assumed for the classical phase-space distribution. As we have 
explained in the previous section, to discuss the correct quantum-mechanical description 
of the problem, the Wigner function must be used instead. This does not invalidate the 
results of |H1 because they are independent of the classical or quantum nature oiW{X, K). 
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Thus counting of coincidences becomes straightforward^^. The number of 
identical events obviously depends on the binning, however, so the procedure 
is ambiguous [ncniiis- To obtain a viable estimate of C{1), we have to 
select the binning in such a way that the result of (jH)) is as close as possible 
to that given by m- 

In jS] the optimal binning was determined and the corresponding relation 
between and C{1) was derived. 

Here we only quote the result in the simplest (but most importantant in 
practice) case when each component of momentum is split into bins of equal 
size (not necessarily the same for each component). Under this condition the 
optimal size of the (3-dimensional) bin is given by 


U 


(^LxLyLz) 


( 16 ) 


With this choice of the binning, the relation between and C{1) is 


C{1) = 


Ebins < > 

Ebins < w{K) 


(17) 


Equations and (HID define the method of estimating the moments of 
Wigner function C{1) from the observed coincidence probabilities One 

sees that, as discussed in detail in jH], the accuracy of the method improves 
for large volume {L^LyL^) of the system. 

It is interesting to observe that the condition (HED for optimal binning 
(i.e. for (fTip to be valid) involves only the product A^AyA^. This may be 
employed to improve the accuracy of the method by selecting small bins in 
the directions where the momentum dependence is strong and large bins in 
the directions whe this dependence is weak^^. This will bring the correction 
factor in (ED closer to 1. 


4 Moments of Wigner function and Renyi en¬ 
tropies 

In this section, using again the general form ED of the Wigner function, 
we discuss the coincidence probabilities C{1) dehned in terms of the density 

detailed description of this procedure was given in [H] and applied in m 
^^E.g., in case of cylindrical symmetry there is no need to split the azimuthal angle. 
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matrix by ©• To this end we observe that, as seen from the density 
matrix for M particles can be expressed as a Fourier transform of Wigner 
function [7j: 


Pm{p]p') = J dXe"^^W{X,K) = 

= (18) 

where q = p — p' and K = {p + p')/2, and where we have explicitely used the 
Gaussian form (fTT|l of G{u). 

When this is introduced into 0 we have 


C{1) = J £‘K,...<pKMa(Ki,-,KM-,l), 


(19) 


where 


0 = / ri Am] + ••• + + ••• + <1m) 


i=l 




yM y [ 

Z_-/7Ti=l = l Z—/a 1 


•'HTna-^^mot {Hma\ ^qlI ^ \ 


and where we have changed the variables from p ^,..., p^^ to 


K = 


X' „(*) ■ „(*) = „(*) _ „(*+!) ■ /yd) = P' 

/ ^ Fm. 5 Hm irm ym ^ o 


d) _i_ T,d+i) 

m 


ItTi 


( 20 ) 

( 21 ) 


with pO+O = p0>. 

One sees from this formula that in the limit of large L only the region 
g ~ 0 contributes signihcantly to the integral. This justihes an expansion of 
v{Ki \ in the exponent. We write 

A'® = A'™ + kS, ( 22 ) 


where are linear combinations of ql^\ This gives up to second order 


!.(A'i‘',.... A'®) = !.(Ai, Km) + v;„lt® + fc®, 


(23) 




where 


^ma 9fno^v(^Ki^ , ^ma.,nl3 ^ma^n/3^(-^ 15 •••5 -^m)* 


(24) 


The indices m, n = 1,..., M denote particles, a, (3 = x,y, z denote directions. 

Introducing (|^ into (PU|) we obtain a gaussian integral which can be 
explicitely evaluated. The result was derived in [TH] and reads 


.... K,, 0 ^ ( Det 


where 


/3A^/2 (L^LyL,)3^(^-i) 

(/-!)! 


1 + 


S = 1 


(25) 


1 


do — 


22 ^ (2s+ !)!(/-2s-1 )! 
and T is the 3M x 3M matrix 


T -J-V — 

-^ma,nl 3 j ^ ma,n 0 j 

Jj(y L/^ 


(26) 


(27) 


Note that the sum over s is hnite, because all vanish for 2s > Z — 1. 

If the eigenvalues of the matrix T are known, the determinant in Eq. 
(ESI) can be explicitely evaluated and one obtains 

iz(Ai,..., Am;Z) - —7777———T 7 —„ f/.—r—r;—r—r— 


(L,4L,)3M(/-i) nEi n„ [(1 + It^aY - (1 + It^aY 

5 Finite volume corrections 

The first thing we observe is that in the limit L —>■ cx) one obtains simply 

(2^)3Ar(/-i)/2 


fZ(Ai,...,i^M) ^3^/2 {L^LyL.f^-^) 

and thus, comparing with da, we have in this limit 

C{1) = 6(1), 


(29) 


( 30 ) 


as expected. 
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One also sees from ((23) that the next order correction depends explicitely 
(through the matrix T) on the shape of the momentum distribution and, 
therefore, it cannot be evaluated in full generality. 

For illustration, we have worked out two examples. 

The simplest case when the particles are uncorrelated and the momentum 
distribution is a Gaussian, 


M 1 

log[2vr/xi] 

m=l a ^ 


(31) 


gives 


and thus 


Tma,n(3 





Det 


1 + Y.^sT^ 

S = 1 



(32) 


(33) 


Consequently, we have 


Ci = 




on 


Oi 


i + E 

S = 1 




(34) 


and one sees that the corrections vanish in the limit —>■ oo. 

As the second example we have studied corrections for / = 3 in a more 
realistic case of an uncorrelated (factorizable) cylindrically symmetric boost- 
invariant distribution with Boltzmann shape in transverse momentum 


M 

m=l 

J K\ + m? , 

v{K) = -+ log E + log[27rT(m + T)e-'"/^] (35) 

where T is a parameter, E is the energy of the particle and m its mass. The 
constant is added to guarantee the correct normalization in the Y interval of 
unit length. Using (123, the matrix Vap and the determinant Det[l -|- T/12], 
necessary to evaluate (|23 for / = 3, can be found [T!)] . 
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Figure 1: The relative difference A 3 = 1 — H{3)/H{3) plotted versus L± for 
several values of the ’’temperature” T. 


In Figure 1, the relative difference 

_ //(3)-g(3) logC(3) 

’’ H{3) logC(3) 

is plotted versus for various values of T. One sees that for T > 150 MeV 
(corresponding to the average transverse momentum larger than 300 MeV), 
and > 3 fm (appropriate for heavy ion collisions) A 3 is indeed very small. 
We thus conclude that the moments of Wigner function reproduce the Renyi 
entropies of a multiparticle system created in high-energy nuclear collisions 
with rather good accuracy. 

One sees also from Fig. 1 that the corrections increase at smaller values 
of L^. At L± 1 fm (appropriate for elementary collisions) they reach 
about 20% for T = 150 MeV and fall quickly at larger T. Thus, although for 
hadron-hadron and lepton-hadron collisions our estimates of Renyi entropy 
clearly require a more precise determination of the size of the system, they 
seem nevertheless to be within reach of the present experiments. 

These two simple examples give only a general idea how to control the 
corrections. Estimates of corrections in more complicated situations are of 
course posssible, e.g., by Monte Carlo simulations. 
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6 Discussion and conclusions 


Several comments are in order. 

(i) Although we have discussed the general case of arbitrary I, it should 
be realized that, in practice, one may at best hope for the determination of 
the lowest order Renyi entropies Z = 2, 3, perhaps also I = 4. This implies 
that an extrapolation to Z = 1, giving the Shannon entropy Q, may require 
an independent input to be reliable ng. 

(ii) One sees from that for Z = 2 ai = 0, i.e., C{2) = 0(2). Thus 
the Renyi entropy H{2) does not suffer from the corrections discussed in this 
paper. 

(iii) As is seen from the discussion in section 4, the difference between 
Renyi entropies and moments of the Wigner function depends primarily on 
the size of the system in configuration space. Also accuracy of the measure¬ 
ments of the moments of Wigner function depends on this size [S], as disussed 
in Section 3. One concludes that information from the HBT measurements 
iniiiHi is important for a successful application of the method. 

(iv) One should keep in mind that we are discussing here the phase-space 
distribution and the Wigner function averaged over time. If the freeze-out 
takes relatively long time, the effective size of the system may be much larger 
than naively expected. On one hand this would improve the accuracy of the 
method. On the other hand it makes the estimate of the volume much more 
difficult, as the standard interpretation of the HBT measurements may be 
not adequate ng. 

(v) It should be emphasized that the estimate of entropy investigated 
in the present paper takes explicitely into account the correlations between 
the particles observed in the experiment. This may be contrasted with the 
estimate obtained in [201 where the entropy of the multiparticle system is 
estimated from the single-particle inclusive distribution, thus ignoring the 
correlations between particles (except those induced by quantum interfer¬ 
ence). It would be very interesting to compare the results obtained from 
these two methods. This may provide an insight into the role of multiparti- 
cle correlations in counting the number of effective degrees of freedom of a 
multiparticle system. 

(vi) The entropies we discuss refer to the particles actually measured in 
experiment in question. Therefore the result does depend on the nature of 
the detector. E.g., the results change when entropy is determined from all 
produced particles instead of only from the charged ones. This should not be 
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surprising: the effective number of degrees of freedom naturally depends on 
the number and nature of the particles considered. Actually, investigations 
of the dependence of entropy on the number of particles may provide some 
interesting hints on the structure of produced final states. 

(vii) Through this paper we have only discussed entropies of the M- 
particle distribution at fixed M. One is often interested in entropies summed 
over all multiplicities. They may be obtained from coincidence probabilities 
referred to all multiplicities, constructed following the formula 

C(l) = i:|P(Af)|'CM(!), (37) 

M 

where P{M) is the multiplicity distribution and where, for the sake of clarity, 
we have added a subspript M to denote the coincidence probability at fixed 
M. One sees that for large I only multiplicities close to the most probable 
one contribute to the sum. 

In conclusion, we have analyzed the relation between the moments of 
Wigner function (jU)) (which can be measured by counting the number of 
identical events jH]) and the coincidence probabilities Q (which define the 
Renyi entropies). It was shown that, for a large class of models, these mo¬ 
ments are identical to the coincidence probabilities in the limit of an infinite 
volume of the system. The finite volume corrections were discussed. They 
were shown to fall as the inverse square of the linear size of the system at 
freeze-out and turn out to be negligible for systems encountered in relativistic 
heavy ion collisions. 
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